Persistent Currents and Dissipation in Narrow Bilayer Quantum Hall Bars 
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Bilayer quantum Hall states support a flow of nearly dissipationless staggered current which can only 
decay through collective channels. We study the dominant finite-temperature dissipation mechanism 
which in narrow bars is driven by thermal nucleation of pseudospin solitons. We find the finite- 
temperature resistivity, predict the resulting staggered current-voltage characteristics, and calculate 
the associated zero-temperature critical staggered current and gate voltage. 
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A 2D electron gas bilayer, subjected to a strong per- 
pendicular magnetic field, can exhibit incompressible 
quantum Hall (QH) states even for filling fractions cor- 
responding to compressible states of noninteracting lay- 
ers The nontrivial, strongly interacting nature of 
these QH states lies in the fact that they survive the limit 
of vanishing interlayer tunneling || . They are stabilized 
by the exchange part of the Coulomb interaction, which, 
in the limit of vanishing single-particle tunneling, sets 
the scale of the gap and leads to macroscopic interlayer 
phase coherence. 

In addition to exhibiting the QHE for a uniform cur- 
rent, these states support persistent currents that are 
counter-propagating in the two layers with J = J top — 
•^bottom- In this Letter, we study a thermally-driven de- 
cay mechanism of J which controls the current-voltage 
characteristics for staggered currents smaller than the 
critical current J c . Because the bilayer system displays a 
quantum Hall gap A in the phase-coherent ground state, 
dissipation via single-particle mechanisms is strongly 
suppressed for fc^T <C A. Therefore, as with supercur- 
rents in superconductors || , the staggered-current decay 
rate is dominated, for a range of parameters, by the col- 
lective mechanism of soliton nucleation. 

A convenient language for describing this strongly cor- 
related quantum-coherent gapped state is in terms of a 
pseudospin unit vector field rh(r) = (m±,m z ) with 
m z = cos 6 = n top — ribottom giving the electron charge- 
density difference between top and bottom layers and 
to_l = sin #(cos</>, sinci) characterizing the relative phase 
4> = 0top — ^bottom of electrons in two layers. The en- 
ergy functional describing long length scale (larger than 
the magnetic length £ — ^/Tic/Be) variations of m(r) is 
given by § 



H 



^\Vm ± \ 2 + ^\Vm z \ 2 + pm 2 z 



(1) 



where electron Coulomb interaction is the origin of the 
effective exchange constants p z s ' that drive the tran- 
sition into the pseudo-ferromagnetic ground state, cor- 
responding to the interlayer phase coherent QH state. 



The electrostatic capacitive energy introduces a hard- 
axis anisotropy, which forces the pseudo-magnetization 
to lie in the _L plane (m z — 0) and thereby reduces the 
full SU(2) pseudospin symmetry to U(l) d^l- A com- 
bination of the external gate voltage V g and the single- 
electron interlayer tunneling t acts as an external pseudo- 
magnetic field h = (t/2iT£ 2 )x + V g z. Because the tunnel- 
ing t can be tuned independently of (3 to be quite small, 
the low-energy physics of this anisotropic QH pseudo- 
ferromagnet, described by the Goldstone mode m, can 
be fully explored experimentally. 

In the limit of vanishing tunneling i, an essentially ex- 
act analytical treatment of narrow (ID limit) QH bars is 
possible and leads to the following results. The bilayer 
QH phase exhibits staggered current-carrying states that 
are metastable and therefore supports staggered persis- 
tent currents for J < J c (V g ), where the critical current 
density is given by 



J c (V g ) = J%(v) 
Z$ jo ^ _ 



(2a) 
(2b) 



where J° = 2y2/3pJ JK is the critical current at zero 
gate voltage, which vanishes in the SU(2) invariant (3 — > 
limit p|, v — V g /2(3 is a reduced measure of the gate volt- 
age V g , and qi{v) = l-w 2 / 3 [( v / 4 + v 2 +2) 1 / 3 - (V4 + v 2 - 
2) 1 / 3 ]. At finite temperature, J c (V g ), plotted in Fig. [I], 
therefore delineates a low resistivity regime, where Ohmic 
dissipation is dominated by slow thermal soliton nucle- 
ation, from a highly resistive state dominated by quasi- 
particle dissipation. 

For J < J c (Vn), the staggered IV characteristics are 
plotted in Fig. and given by 



V(V g ,J) 



hujoL 

e£ 



sinh (7T7'T /2T) e 



-(U B +ivj/2)T /T 



(3) 



where luq = \j2ixP-tf3 jh is the microscopic attempt fre- 
quency, j = J/ J® is the reduced (dimensionless) current 
density, fc^To = hJ®L y , with L y the narrow sample di- 
mension <C L x = L, and Eb = ksToUB is the saddle- 
point energy barrier separating two different current- 
carrying states. The barrier is plotted in the inset of 
Fig. and is explicitly given by 
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FIG. 1. Main Plot: Critical staggered current j c as a func- 
tion of gate voltage v. Upper Inset: Sample geometry consid- 
ered in this work. Lower Inset: Pseudospin configuration of 
the 27r-line soliton whose nucleation leads to staggered current 
decay. 



FIG. 2. Main Plot: IV curve for J < J c (V g ), pi = p^, 
V g = and To/T = 6. Inset: Energy barrier for various gate 
voltages v = V g /2/3, with the barrier vanishing for v 1. 
Qualitatively similar curves are obtained for the full range 
< p z s /pi < 1. 



U B (Vg,J) 



pg- 1 — m 2 
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+ 2v (^Jl-m 2 , - \J\ - m 2 ^ 



1/2 

, (4) 



where the limits of integration are to 2 , = j/q and mf = 
q 2 [m 2 l + 2fh 2 (l + y/l + ml/fh 2 )], with m 2 = (l-mg)(l- 
q ), and the dimensionless wavevector q is defined implic- 
itly through the current j = q[l — v 2 /(l — q 2 ) 2 ]. The 
analytic expression for the barrier simplifies consider- 
ably when one of its arguments vanishes. We find that 
U B (V g ,J) depends weakly on e = 1 — (pf/p^) and for 
e — ^ 1 is given by 



U B (0,J) 
U B (V g ,0) 
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(5a) 
(5b) 



and U B (0,0) = \-JT~£~ 1 ±j t (| -sin" 1 (l-2e)). For 
such narrow Hall bars the staggered linear resistivity p is 
always finite at finite temperature and is given by 



P(Vg) 



V j^o h 2 u>oL 



LJe 



2k B Te 2 £ 



v e -u B (v g ,o)T /T_ 



(6) 



From exact diagonalization studies g] at d/i = 1/2, we 
have = 5.1 X lCT 3 e 2 /^ 3 e and pj- = 1.5 x lCT 2 e 2 /£e. 
Taking in addition £ = 20 nm, T = 1 K, e = 13.2, t 
0.1 mrY. = 0, and e = 0, we obtain 



where iV = L y /£, Setting N = 5, 10, 15 respectively 
gives p(0) = 245 ft, 0.217, and 115 /ifi. 

For a realistic system, there is a limited range of 
validity of the above results, with other effects domi- 
nating outside of this range. The constraint of quasi- 
equilibrium, which implies low decay rate, together with 
the requirement that the thermal collective dissipation 
mechanism dominates over single-particle current decay 
requires k B T < E B < A. At the same time, however, 
T must be sufficiently high so that thermal nucleation 
dominates over quantum tunneling of phase slips, j^] Fur- 
thermore, in order for the bulk nucleation rate to be 
experimentally observable, it is necessary that it dom- 
inates over phase slips nucleated at surfaces, contacts, 
and sample inhomogeneities. Since the bulk nucleation 
rate scales with the Hall bar length L, we expect that 
the bulk mechanism dominates over surface nucleation 
for L ^> £. Also, for the staggered current decay rate to 
be dominated by the ID line solitons studied here, the 
saddle-point energy barrier given in Eq. (Q) must be lower 
than barriers for the competing mechanism of ± vortex 
pair nucleation. For a sufficiently wide Hall bar, the lat- 
ter scenario will dominate, with the crossover occurring 
for L y « £, = ^/pjjp. 

We now present the highlights of calculations that lead 
to these results. Although quite distinct in detail, the 
spirit of our analysis follows the classic work of Langer 
and Ambegaokar j|. 

The Euler-Lagrange (EL) equations for Eq. ([!]) admit 
the uniform current-carrying solutions 



p(o) « io 5 iv(o.2) JV n 



(7) 



4>{k) = kx, 
m 2 ± (v,q) = l-v 2 /(l-q 2 ) 2 , 



(8a) 
(8b) 
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FIG. 3. Main Plot: Effective potential, Eq. (10b), for the 
mechanical analogy of the EL equations plotted for v = 0.5 
and j — 0.1. The physically stable configuration is the me- 
chanically unstable point Inset: Schematic drawing of 
the saddle-point (bounce) solution — the dominant contribu- 
tion to the transition probability. 



where q 2 = k 2 pjr/2(3 is the dimensionless wave vector 
and we have taken our sample to lie in the x-y plane 
with dimensions L x = L L y . Equation ( |Sb| ) is valid in 
the region q 2 < 1 — v and v = V g /2(3 < 1. (We consider 
V g to be non-negative.) The staggered current for this 
solution is J = 2p^m 2 1 k/h, or, equivalently, j — ra^q. 

For nonuniform solutions, the EL equations can be 
combined into a single equation which, after some ma- 
nipulation, can be written as 



X -M (d x m ± ) 2 + V cS = E eS , 

for some constant E c g, with 
M(m ± ) =p^+p>l/(l-mi), 
V cB (m x ) = P ( ~^~T ~ i 1 ~ m D + 2u \A 



(9) 



(10a) 
(10b) 



In the usual mechanical analogy, Eq. (^) represents the 
energy E c g of a particle at "position" m± and "time" 
x, moving in a potential V c g(m±) and with a space- 
dependent mass M(m±). This potential is plotted in 
Fig. | for v = 0.5 (j c re 0.28) and j = 0.1. 

The "conservation of energy," Eq. , immediately im- 
plies the existence of two extended solutions || . A uni- 
form current-carrying solution is given by Eq. cor- 
responding to the particle forever remaining on top of 
the hill at m 1 ^ with angular velocity k. For this so- 
lution, winds azimuthally at a constant rate as a 
function of x with a constant amplitude deviating 
from the equatorial plane (mj_ = 1) with increasing V g . 
For fixed winding k, this uniform current-carrying solu- 
tion is a local minimum of H and is hence metastable. 



Since the energy is lowered with decreasing k, we ex- 
pect at finite temperature the system will thermally ac- 
tivate down to the k = zero-current ground state via 
successive thermal transitions k — > k — 2it/L. To cal- 
culate the rate of such transitions we need to compute 
the saddle-point barrier separating two "neighboring" 
current-carrying metastable states. Brief reflection on 
the mechanical analogy shows that the second solution is 
in fact a saddle-point solution corresponding to the par- 
ticle starting at "position" and "time" x — — oo, 
spiraling down hill to rrS^ at x — while conserving 
angular momentum but increasing its angular velocity 
k(m±) oc l/m\, and finally bouncing back out to 
for x — +oo. The resulting "energy conservation" equa- 
tion allows an exact determination of the saddle-point 
solution, written in terms of a ID integral, 



1 



m± (x) 



dm ± yjM( mi _)/[E c ff— V e s(m±)], (11) 



for an appropriate value of E e g, and is shown in the inset 
of Fig. | 

The saddle-point solution is the nucleation site for 
the eventual singularity where m± — > (m z — ► ±1), 
and whereby the system can slip a loop reducing the 
total phase winding Aip by 2ir, and therefore leading 
to staggered-current decay. Contrary to what is of- 
ten tacitly assumed, the barrier for defect nucleation 
is not determined by the energy of the defect — which 
may be either singular, as it is in superconductors, or 
non-singular, as it is in the present case of the full 
SU(2) space — although we expect them to be close in en- 
ergy. Instead, the barrier is controlled by a non-singular 
saddle-point field configuration (inset of Fig. ||), which 
is in the same topological sector as the current-carrying 
metastable minimum. It is also important to note that 
for sufficiently small /?, the energy for such phase slips 
is much smaller than the quantum Hall gap A and that 
therefore the entire system remains in the fully gapped 
quantum Hall state throughout the phase-slip process. 
This is to be contrasted with phase slips inside a super- 
conductor, where the order parameter, and therefore the 
bulk gap, are suppressed within the vortex core and one 
in principle has to take into account the core-confined 
low energy (normal) quasi-particle degrees of freedom. 

The energy barrier is defined as the difference in en- 
ergy between the saddle-point solution, Eq. ( |Tl| ) and the 
uniform current-carrying solution, Eq. (^). By exploiting 
the mechanical analogy, and in particular the conserva- 
tion law, Eq. (|^), we obtain the expression quoted in 
Eq. |). 

In the steady-state regime, the staggered voltage 
V(V g ,J) is proportional to the net rate of phase 
slips, c>tA0/27r, which itself is the difference be- 
tween the rate of current-decreasing transitions (= 
uiq exp[— EB/kBT]) and the rate of current-increasing 



3 



transitions (= uq exp[— (Eb + irhJ L y ) /ksT]). In a sam- 
ple of length L, there are approximately L/£ possible nu- 
cleation sites. These considerations lead directly to the 
expression in Eq. (||). 

In contrast to the energy barrier Eb , which is a static 
quantity, we must consider pseudospin dynamics in or- 
der to compute the attempt frequency ujq appearing in 
Eq. (|). Within the microscopic dynamical model valid 
at low temperatures, T <C t, for finite interlayer tun- 
neling, ujq is given by the ratio of the curvature of the 
metastable well (~ t/2ir£ 2 ) and the dynamical mass term 
(~ ^ 2 /4tt 2 £ 4 /3), leading to hu) = yj2n£ 2 t{3. In contrast, 
for T t, we expect classical Langevin dynamics ||, 
characterized by a kinetic "drag" coefficient 7, and which 
in the simplest estimate gives ujq = /3/7- 

When the transverse (narrow) dimension L y of the 
sample becomes sufficiently large, the energy barrier 
Eb ~ (pjrPY^Ly for nucleating a line-soliton defect be- 
comes comparable to the energy of nucleating a ± vortex 
pair — a competing mechanism for inducing phase-slips. 
Up to weak logarithmic corrections, the energy of such a 
vortex pair is E v w pj: — (3S y 2 , where £ = (pj'/P) 1 ^ 2 is the 
core size. Vortex nucleation should therefore dominate 
the ID soliton nucleation considered here for L y > £, 
which can be tuned independently of the QH gap A. 
This is in contrast with superconductors, where the cor- 
responding ID to 2D crossover scale is the Ginzburg- 
Landau correlation length, controlled by the supercon- 
ducting gap. 

Although it is difficult to extend our exact ID analysis 
to the 2D limit, we can estimate the staggered current 
decay rate using simple scaling arguments. In 2D, the 
phase slip rate is controlled by ± vortex pair nucleation, 
analogously to superfluids and superconductors. How- 
ever, in contrast to those more familiar systems, here 
vortices (half-skyrmion, i.e., merons) carry ±1/2 elec- 
tromagnetic charge in addition to their U(l) topological 
charge, and there are therefore four elementary vortex 
defects (+1/2, L), (—1/2, L), (+1/2, R), and (-1/2, R), 
with L, R respectively corresponding to ±27r circulation 
of <j). Correspondingly, within the interlayer phase (</>) 
coherent state, the L and R vortices are bound into 2 
types of topologically neutral pairs: (i) electromagnet- 
ically neutral pairs [(+1/2, L)-{— 1/2, R)} and (ii) elec- 
tromagnetically charged pairs [(+1/2, L)-(+l/2, R)] 0. 
Nevertheless, we do not expect the Coulomb interaction, 
which is subdominant to the topological-charge confin- 
ing potential, to play a role in staggered-current-induced 
vortex ionization processes. Hence, in the limit of van- 
ishing t, standard nucleation analysis for dissipation in 
superconducting films j|] can be easily extended to our 
system. It predicts a highly nonlinear power-law stag- 
gered IV, E ~ J Q , with a(T) > 3 in the interlayer co- 
herent state, a result that contrasts strikingly with the 
linear staggered IV found in the ID limit. 



Nonvanishing interlayer tunneling, t explicitly breaks 
C/(l) symmetry and leads to nonuniform staggered 
current-carrying states composed of a lattice of solitons 
of width S = £(2np±/t) 1 / 2 and density n akin to a peri- 
odic array of Rayleigh-Benard current rolls. |9|. 

In the dense soliton limit, n5 ^> 1, (relevant for small 
t and large J) the current J(x) is nearly uniform and 
our ID results directly apply. In the dilute limit, nS ^> 1 
(which is always reached for sufficiently low J > J c \ ~ \fi 
pO|,p[), for £ <C 5, phase slips are confined to a single 
soliton, inside which d x <f> « 2it /5 is uniform and the tun- 
neling energy is on average zero. Hence, our ID, t = 
analysis again applies with the J-independent wavevec- 
tor k cS = 2ir/5. Furthermore, scaling analysis |J sug- 
gests that in the opposite limit £ > <5, this solution is still 
valid, but with the nucleation width set by 5, rather than 

In the 2D limit, the staggered current decays by ion- 
ization of L-R vortex pairs, whose energy for t > grows 
linearly with the separation R as (pstY^R/i. This there- 
fore suggests the existence of a true staggered critical 
current J c = (e/h)y/p s t/£ 2 , with E( J < J c ) = even at 
finite T (up to corrections that vanish in the thermody- 
namic limit), and for J > J c , E(J) « | J- J c | q(t) . 
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